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Q_i Abstract 

^ Measurements of b hadron production ratios in proton-proton collisions at a centre- 

' of- mass energy of 7 TeV with an integrated luminosity of 3 pb~^ are presented. We 

study the ratios of strange B meson to light B meson production fs/{fu + fd) and 
baryon to light B meson production fA^/ifu + fd) as a function of the charmed 
hadron-muon pair transverse momentum and the b hadron pseudorapidity rj, 
^ for pt between and 14 GeV and rj between 2 and 5. We find that fs/ (/« -|- fd) 

CN is consistent with being independent of px and r/, and we determine fs/{fu + fd) 

where the first error is statistical and the second systematic. 
^ The corresponding ratio fA^/{fu + fd) is found to be dependent upon the transverse 

^ momentum of the charmed hadron-muon pair, fK^/{fu + fd) = (0.404ib0.017(stat)± 

0.027(syst)±0.105(Br)) x [l-(0.031±0.004(stat)±0.003(syst)) xpT(GeV)], where Br 
refiects an absolute scale uncertainty due to the poorly known branching fraction 
^ B{A^ — )• pK~TT~^). We extract the ratio of strange B meson to light neutral B 

meson production fs/fd by averaging the result reported here with two previous 
measurements derived from the relative abundances of B^ — )• Df-K^ to B^ — J- D^K~ 
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and 5° ^ D+-K-. We obtain fs/fd = 0.267+° °^^ 
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1 Introduction 

The fragmentation process, in which a primary b quark forms either a bq meson or a 
bqiq2 baryon, cannot be rehably predicted because it is driven by strong dynamics in the 
non-perturbative regime. Thus fragmentation functions for the various hadron species 
must be determined experimentally. The LHCb experiment at the LHC explores a unique 
kinematic region: it detects b hadrons produced in a cone centered around the beam axis 
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covering a region of pseudorapidity rj, defined in terms of the polar angle 6 with respect to 
the beam direction as — ln(tan^/2), ranging approximately between 2 and 5. Knowledge 
of the fragmentation functions allows us to relate theoretical predictions of the bb quark 
production cross-section, derived from perturbative QCD, to the observed hadrons. In 
addition, since many absolute branching fractions of B~ and decays have been well 
measured at e+e~ coUiders yy, it suffices to measure the ratio of -B° production to either 
or 5° production to perform precise absolute B^ branching fraction measurements. In 
this paper we describe measurements of two ratios of fragmentation functions: fs/ {fu + fd) 
and /a^/ {fu + fd), where fq = B{b — )■ Bq) and /a^ = B{b — )■ A;,). The inclusion of charged 
conjugate modes is implied throughout the paper, and we measure the average production 
ratios. 

Previous measurements of these fractions have been made at LEP [2j and at CDF [H] . 
More recently, LHCb measured the ratio fs/fd using the decay modes B^ — )■ D^n^ , B^ — t- 
D^K~ , and 5° — )■ Df7T~ [1] and theoretical input from QCD factorization [5l|6]. Here we 
measure this ratio using semileptonic decays without any significant model dependence. 
A commonly adopted assumption is that the fractions of these different species should be 
the same in high energy b jets originating from decays and high b jets originating 
from pp collisions at the Tevatron or pp collisions at LHC, based on the notion that 
hadronization is a non-perturbative process occurring at the scale of Aqcd- Nonetheless, 
the results from different experiments are discrepant in the case of the b baryon fraction [2] . 

The measurements reported in this paper are performed using the LHCb detector [7J , 
a forward spectrometer designed to study production and decays of hadrons containing 
6 or c quarks. LHCb includes a vertex detector (VELO), providing precise locations 
of primary pp interaction vertices, and of detached vertices of long lived hadrons. The 
momenta of charged particles are determined using information from the VELO together 
with the rest of the tracking system, composed of a large area silicon tracker located 
before a 4 Tm dipole magnet, and a combination of silicon strip and straw drift chamber 
detectors located after the magnet. Two Ring Imaging Cherenkov (RICH) detectors are 
used for charged hadron identification. Photon detection and electron identification are 
implemented through an electromagnetic calorimeter followed by a hadron calorimeter. 
A system of alternating layers of iron and chambers provides muon identification. The 
two calorimeters and the muon system provide the energy and momentum information 
to implement a first level (LO) hardware trigger. An additional trigger level is software 
based, and its algorithms are tuned to the experiment's operating condition. 

In this analysis we use a data sample of 3 pb~^ collected from 7 TeV centre-of-mass 
energy pp collisions at the LHC during 2010. The trigger selects events where a single 
muon is detected without biasing the impact parameter distribution of the decay products 
of the b hadron, nor any kinematic variable relevant to semileptonic decays. These features 
reduce the systematic uncertainty in the efficiency. Our goal is to measure two specific 
production ratios: that of -B^ relative to the sum of B~ and B^, and that of A'^, relative 
to the sum of B~ and B^. The sum of the B^, B~ , B^ and fractions does not equal 
one, as there is other b production, namely a very small rate for B~ mesons, bottomonia, 
and other b baryons that do not decay strongly into such as the Eh- We measure 
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Table 1: Charmed hadron decay modes and branching fractions. 



Particle 


Final state 


Branching fraction (%) 






3.89±0.05 [U 






9.i4±o.2o ng 


Dt 




5.50±0.27 [20] 






5.0±1.3 pp 



relative fractions by studying the final states D^fi~uX, D^fi~VX, n~VX, yl+/i~z7X, 
D^K~^ fi~VX , and D^pfi~uX. We do not attempt to separate fu and fa, but we measure 
the sum of D° and channels and correct for cross-feeds from and decays. We 
assume near equality of the semileptonic decay width of all b hadrons, as discussed below. 
Charmed hadrons are reconstructed through the modes listed in Table [T| together with 
their branching fractions. We use all Df — )■ K~K~^tt~^ decays rather than a combination 

of the resonant ^vr"*" and K contributions, because these decays cannot be cleanly 
isolated due to interference effects of different amplitudes. 

Each of these different charmed hadron plus muon final states can be populated by a 
combination of initial h hadron states. mesons decay semileptonically into a mixture 
of and mesons, while B~ mesons decay predominantly into mesons with a 
smaller admixture of mesons. Both include a tiny component of D^K meson pairs. 

mesons decay predominantly into mesons, but can also decay into D^K^ and D~^Kg 
mesons; this is expected if the B^ decays into a D** state that is heavy enough to decay 
into a DK pair. In this paper we measure this contribution using D^K~^X fi~V events. 
Finally, baryons decay mostly into A'^ final states. We determine other contributions 
using D^pX^~V events. We ignore the contributions of 6 — > m decays that comprise 
approximately 1% of semileptonic h hadron decays [8J, and constitute a roughly equal 
portion of each h species in any case. 

The corrected yields for _B° or B~ decaying into D^ii'VX or D^fi^VX, ricovi, can be 
expressed in terms of the measured yields, n, as 



i3(D0 ^ K-n+)e{B 



X 



(1) 



e(ylo ^ D^p) 



where we use the shorthand n{Dfi) = n{DX fi~V). An analogous abbreviation e is used for 
the total trigger and detection efficiencies. For example, the ratio e(i?° — )■ D°)/e(i?° — )■ 
D^K^) gives the relative efficiency to reconstruct a charged K in semi-muonic 5° decays 
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producing a meson. Similarly 



e{B D+) [B{D+ K-7r+7T+) 

n{D^K+fi-) e(50 ^ D+) 
B{DO K-TT+) e(50 ^ DOA'+) 

B{D^ K-71+) e{Ab DOp)\ ' ^ ^ 

Both the D^Xn'T/ and the D^Xfi^u final states contain small components of cross-feed 
from i?g decays to D^K^ X ix'V and to D~^K^Xfi~u. These components are accounted for 
by the two decays — t- DfiXfi^V and — t- D*^X/i~z7 as reported in a recent LHCb 
publication [9]. The third terms in Eqs. [l] and [2] are due to a similar small cross-feed from 
yl° decays. 

The number of resulting in DfXfi'u in the final state is given by 



B{D+ K+K-n- 



N{B^ + B-)B{B DjKn)e{B DfKfx 



where the last term subtracts yields of D^KXfi v final states originating from B^ or B 
semileptonic decays, and N{B^ + B~) indicates the total number of B^ and B~ produced. 
We derive this correction using the branching fraction B{B — )■ Di*^~^ K nu) = (6.1 ± 1.2) x 
10~^ [T^ measured by the BaBar experiment. In addition, B^ decays semileptonically 
into DRXfi'V, and thus we need to add to Eq. Is] 



neo„(5° ^ DKfx) = 2 niD'^K+ji) 

where, using isospin symmetry, the factor of 2 accounts for B^ DK^Xfi'u semileptonic 
decays. 

The equation for the ratio fs/{fu + fd) is 



fu + fd '^corr 

{B DV) + '^corr 



(5) 



where B^ ^ represents B^ semileptonic decays to a final charmed hadron, given 
by the sum of the contributions shown in Eqs. [3] and [4| and the symbols r^- indicate 
the Bi hadron lifetimes, that are all well measured pj. We use the average B^ lifetime, 
1.472±0.025 ps [1]. This equation assumes equality of the semileptonic widths of all the b 
meson species. This is a reliable assumption, as corrections in HQET arise only to order 
I /ml and the SU(3) breaking correction is quite small, of the order of 1% [HI [T21 [IS]- 
The corrected yield is derived in an analogous manner. We determine 

(AO.n\- njAtfj-) n{D°p^i~) 
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where D represents a generic charmed hadron, and extract the A° fraction using 

Ik ^corr(Aj; ^ D/i) Tb- + ^gO q c\ 

Again, we assume near equahty of the semileptonic widths of different b hadrons, but 
we apply a small adjustment ^ = 4±2%, to account for the chromomagnetic correction, 
affecting 6-flavoured mesons but not b baryons [HI |12l |13]. The uncertainty is evaluated 
with very conservative assumptions for all the parameters of the heavy quark expansion. 



2 Analysis method 

To isolate a sample of b flavoured hadrons with low backgrounds, we match charmed 
hadron candidates with tracks identified as muons. Right-sign (RS) combinations have 
the sign of the charge of the muon being the same as the charge of the kaon in D^, 
D~^, or /l+ decays, or the opposite charge of the pion in decays, while wrong-sign 
(WS) combinations comprise combinations with opposite charge correlations. WS events 
are useful to estimate certain backgrounds. This analysis follows our previous inves- 
tigation of 6 —7- D^X^~V [H]. We consider events where a well-identified muon with 
momentum greater than 3 GeV and transverse momentum greater than 1.2 GeV is found. 
Charmed hadron candidates are formed from hadrons with momenta greater than 2 GeV 
and transverse momenta greater than 0.3 GeV, and we require that the average transverse 
momentum of the hadrons forming the candidate be greater than 0.7 GeV. Kaons, pions, 
and protons are identified using the RICH system. The impact parameter (IP), defined 
as the minimum distance of approach of the track with respect to the primary vertex, 
is used to select tracks coming from charm decays. We require that the x^? formed by 
using the hypothesis that each track's IP is equal to 0, is greater than 9. Moreover, the 
selected tracks must be consistent with coming from a common vertex: the per number 
of degrees of freedom of the vertex fit must be smaller than 6. In order to ensure that the 
charm vertex is distinct from the primary pp interaction vertex, we require that the x^i 
based on the hypothesis that the decay flight distance from the primary vertex is zero, is 
greater than 100. 

Charmed hadrons and muons are combined to form a partially reconstructed b hadron 
by requiring that they come from a common vertex, and that the cosine of the angle 
between the momentum of the charmed hadron and muon pair and the line from the 
vertex to the primary vertex be greater than 0.999. As the charmed hadron is a decay 
product of the b hadron, we require that the difference in z component of the decay vertex 
of the charmed hadron candidate and that of the beauty candidate be greater than 0. We 
explicitly require that the rj of the b hadron candidate be between 2 and 5. We measure 
r] using the line defined by connecting the primary event vertex and the vertex formed by 
the D and the /x. Finally, the invariant mass of the charmed hadron and muon system 
must be between 3 and 5 GeV for D^ix~ and D~^fj,~ candidates, between 3.1 and 5.1 GeV 
for D^fi~ candidates, and between 3.3 and 5.3 GeV for A'^fi~ candidates. 
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We perform our analysis in a grid of 3 77 and 5 bins, covering the range 2 < rj < 5 
and pt < 14 GeV. The b hadron signal is separated from various sources of background 
by studying the two dimensional distribution of charmed hadron candidate invariant mass 
and ln(IP/mm). This approach allows us to determine the background coming from false 
combinations under the charmed hadron signal mass peak directly. The study of the 
ln(IP/mm) distribution allows the separation of prompt charm decay candidates from 
charmed hadron daughters of b hadrons [H]. We refer to these samples as Prompt and 
Dfb respectively. 

2.1 Signal extraction 

We describe the method used to extract the charmed hadron-/i signal by using the 
D^Xfi~17 final state as an example; the same procedure is applied to the final states 
D^Xfi~iJ, DfXfi~lJ, and A'^Xfi'u. We perform unbinned extended maximum likelihood 
fits to the two-dimensional distributions in K~tt~^ invariant mass over a region extending 
ibSO MeV from the mass peak, and ln(IP/mm). The parameters of the IP distribution 
of the Prompt sample are found by examining directly produced charm [13] whereas a 
shape derived from simulation is used for the Dfb component. 

An example fit for D^ijTVX, using the whole px and rj range, is shown in Fig. [l] 
The fitted yields for RS are 27666±187 Dfb, 695±43 Prompt, and 1492±30 false 
combinations, inferred from the fitted yields in the sideband mass regions, spanning the 
intervals between 35 and 75 MeV from the signal peak on both sides. For WS we find 
362±39 Dfb, 187±18 Prompt, and 1134±19 false D° combinations. The RS yield includes 
a background of around 0.5% from incorrectly identified /i candidates. As this paper 
focuses on ratios of yields, we do not subtract this component. Figure [2] shows the 
corresponding fits for the D^X^'V final state. The fitted yields consist of 9257±110 Dfb 
events, 362±34 Prompt, and 1150±22 false combinations. For WS we find 77±22 
Dfb, 139±14 Prompt and 307±10 false combinations. 

The analysis for the D^Xfi~u mode follows in the same manner. Here, however, we 
are concerned about the refiection from — )■ pK~7T~^ where the proton is taken to be 
a kaon, since we do not impose an explicit proton veto. Using such a veto would lose 
30% of the signal and also introduce a systematic error. We choose to model separately 
this particular background. We add a probability density function (PDF) determined 
from simulation to model this, and the level is allowed to fioat within the estimated error 
on the size of the background. The small peak near 2010 MeV in Fig. ^h) is due to 
D*~^ — )■ 7j'~^D^,D^ — )■ K^K~. We explicitly include this term in the fit, assuming the 
shape to be the same as for the signal, and we obtain 4±1 events in the RS signal 
region and no events in the WS signal region. The measured yields in the RS sample are 
2192±64 Dfb, 63±16 Prompt, 985±145 false Df background, and 387±132 A+ refiection 
background. The corresponding yields in the WS sample are 13±19, 20±7, 499±16, and 
3±3 respectively. Figure [3] shows the fit results. 

The last final state considered is A'^X^'V. Figure |4] shows the data and fit components 
to the ln(IP/mm) and pK^TC^ invariant mass combinations for events with 2 < rj < 5. 
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Figure 1: The logarithm of the IP distributions for (a) RS and (c) WS candidate 
combinations with a muon. The dotted curves show the false background, the small 
red-solid curves the Prompt yields, the dashed curves the Dfb signal, and the larger green- 
solid curves the total yields. The invariant K'tt'^ mass spectra for (b) RS combinations 
and (d) WS combinations are also shown. 

This fit gives 3028±112 RS Dfb events, 43±17 RS Prompt events, 589±27 RS false A+ 
combinations, 9±16 WS Dfb events, 0.5±4 WS Prompt events, and 177±10 WS false A'^ 
combinations. 

The may also decay into D^pX fi^u. We search for these decays by requiring the 
presence of a track well identified as a proton and detached from any primary vertex. 
The resulting D^p invariant mass distribution is shown in Fig. [sj We also show the 
combinations that cannot arise from decay, namely those with D^p combinations. 
There is a clear excess of RS over WS combinations especially near threshold. Fits to 
the K~'n'~^ invariant mass in the [m{K^7i^p) — m{K^T{^) + m{D^)-pDQ\ region shown in 
Fig. [sj^a) give 154±13 RS events and 55±8 WS events. In this case, we use the WS yield 
for background subtraction, scaled by the RS/WS background ratio determined with a 
MC simulation including {B~ + ^ D^X^~v) and generic hb events. This ratio is 
found to be 1.4±0.2. Thus, the net signal is 76±17±11, where the last error reflects the 
uncertainty in the ratio between RS and WS background. 
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Figure 2: The logarithm of the IP distributions for (a) RS and (c) WS candidate 
combinations with a muon. The grey-dotted curves show the false background, the 
small red-solid curves the Prompt yields, the blue-dashed curves the Dfb signal, and the 
larger green-solid curves the total yields. The invariant K~tt^tt^ mass spectra for (b) RS 
combinations and (d) WS combinations are also shown. 
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Figure 3: The logarithm of the IP distributions for (a) RS and (c) WS candidate 
combinations with a muon. The grey-dotted curves show the false Df background, the 
small red-solid curves the Prompt yields, the blue-dashed curves the Dfb signal, the 
purple dash-dotted curves represent the background originating from A'^ reflection, and 
the larger green-solid curves the total yields. The invariant K^K^tt^ mass spectra for 
RS combinations (b) and WS combinations (d) are also shown. 



9 




ln(IP/mm) m(pK"7iT (MeV) 

Figure 4: The logarithm of the IP distributions for (a) RS and (c) WS candidate 
combinations with a muon. The grey-dotted curves show the false background, the 
small red-solid curves the Prompt yields, the blue-dashed curves the Dfb signal, and the 
larger green-solid curves the total yields. The invariant pK~7r~^ mass spectra for RS 
combinations (b) and WS combinations (d) are also shown. 
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Figure 5: (a) Invariant mass of D^p candidates that vertex with each other and together 
with a RS muon (black closed points) and for a p (red open points) instead of a p; (b) fit 
to invariant mass for RS events with the invariant mass of D'^p candidate in the signal 
mass difference window; (c) fit to invariant mass for WS events with the invariant 
mass of D^p candidate in the signal mass difference window. 
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Figure 6: Projections of the two-dimensional fit to the and m{D^fi) distributions 
of semileptonic decays including a Df meson. The D*/Ds ratio has been fixed to the 
measured D*/D ratio in light B decays (2.42±0.10), and the background contribution is 
obtained using the sidebands in the K'^K^tt^ mass spectrum. The different components 
are stacked: the background is represented by a black dot-dashed line, by a red dashed 
hne, D*^ by a blue dash-double dotted line and D**^ by a green dash-dotted line. 

2.2 Background studies 

Apart from false D combinations, separated from the signal by the two-dimensional fit 
described above, there are also physical background sources that affect the RS Dfb sam- 
ples, and originate from bb events, which are studied with a MC simulation. In the meson 
case, the background mainly comes from b — )■ DDX with one of the D mesons decay- 
ing semi-muonically, and from combinations of tracks from the pp — >■ bbX events, where 
one b hadron decays into a D meson and the other b hadron decays semi-muonically. 
The background fractions are (1.9±0.3)% for D^Xfi~V, (2.5±0.6)% for D~^Xfi~lJ, and 
(5.1±1.7)% for D^XijTV. The main background component for semileptonic decays 
is decaying into D~ A'^, and the D~ decaying semi-muonically. Overall, we find a very 
small background rate of (1.0±0.2)%, where the error reflects only the statistical uncer- 
tainty in the simulation. We correct the candidate b hadron yields in the signal region 
with the predicted background fractions. A conservative 3% systematic uncertainty in 
the background subtraction is assigned to reflect modelling uncertainties. 

2.3 Monte Carlo simulation and efficiency determination 

In order to estimate the detection efficiency, we need some knowledge of the different 
final states which contribute to the Cabibbo favoured semileptonic width, as some of the 
selection criteria affect final states with distinct masses and quantum numbers differently. 
Although much is known about the and B~ semileptonic decays, information on the 
corresponding 5° and A^ semileptonic decays is rather sparse. In particular, the hadronic 
composition of the final states in 5° decays is poorly known [H], and only a study from 



12 



CDF provides some constraints on the branching ratios of final states dominant in the 
corresponding yl° decays fTE]. 

In the case of the -B^ — )■ semileptonic decays, we assume that the final states are 
D+, D*+, L)*o(2317)+, L>,i(2460)+, and D,i(2536)+. States above DK threshold decay 
predominantly into D^*^K final states. We model the decays to the final states D^/i~z7 and 
D*~^fi~u with HQET form factors using normalization coefficients derived from studies 
of the corresponding and B~ semileptonic decays [T|, while we use the ISGW2 form 
factor model p3] to describe final states including higher mass resonances. 

In order to determine the ratio between the different hadron species in the fi- 
nal state, we use the measured kinematic distributions of the quasi-exclusive process 
5° — )■ D^^~VX. To reconstruct the squared invariant mass of the pair (g^), we 
exploit the measured direction of the h hadron momentum, which, together with energy 
and momentum conservation, assuming no missing particles other than the neutrino, al- 
low the reconstruction of the p 4-vector, up to a two-fold ambiguity, due to its unknown 
orientation with respect to the B flight path in its rest frame. We choose the solution cor- 
responding to the lowest h hadron momentum. This method works well when there are no 
missing particles, or when the missing particles are soft, as in the case when the charmed 
system is a D* meson. We then perform a two-dimensional flt to the versus m{fiDf) 
distribution. Figure |6] shows stacked histograms of the Df, D*^, and D**^ components. 
In the flt we constrain the ratio i3(-B° — > D*'^ ^~V)/B{B'^ — )■ ii~V) to be equal to the 
average D*fi~u/Dn~u ratio in semileptonic B^ and B~ decays (2.42±0.10) [1]. This con- 
straint reduces the uncertainty of one D** fraction. We have also performed flts removing 
this assumption, and the variation between the different components is used to assess the 
modelling systematic uncertainty. 

A similar procedure is applied to the A^^~ sample and the results are shown in Fig. [7j 
In this case we consider three flnal states, A^fi^V, Ac{2595)~^ fi~V, and A c (2625)~^ fi~T', with 
form factors from the model of Ref. [17]. We constrain the two highest mass hadrons to 
be produced in the ratio predicted by this theory. 

The measured pion, kaon and proton identiflcation efficiencies are determined using 
K^, D*^, and A° calibration samples where p, K, and tt are selected without utilizing 
the particle identification criteria. The efficiency is obtained by fitting simultaneously the 
invariant mass distributions of events either passing or failing the identification require- 
ments. Values are obtained in bins of the particle t] and p^, and these efficiency matrices 
are applied to the MC simulation. Alternatively, the particle identification efficiency can 
be determined by using the measured efficiencies and combining them with weights pro- 
portional to the fraction of particle types with a given r] and p^ for each /i charmed hadron 
pair 7] and p^ bin. The overall efficiencies obtained with these two methods are consistent. 
An example of the resulting particle identification efficiency as a function of the rj and p^ 
of the A^n~ pair is shown in Fig. |8} 

As the functional forms of the fragmentation ratios in terms of p^ and rj are not known, 
we determine the efficiencies for the final states studied as a function of pt and t] within 
the LHCb acceptance. Figure [9] shows the results. 
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Figure 7: Projections of the two-dimensional fit to the and m{A^ix^) distributions 
of semileptonic decays including a baryon. The different components are stacked: 
the dotted line represents the combinatoric background, the bigger dashed line (red) 
represents the n~V component, the smaller dashed line (blue) the /lc(2595)"'", and the 
solid line represents the ylc(2625)^ component. The ylc(2595)"''//lc(2625)"'' ratio is fixed 
to its predicted value, as described in the text. 
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Figure 8: Measured proton identification efficiency as a function of the A'^fi~ for 
2<?7<3, 3<r7<4, 4<r7<5 respectively, and for the selection criteria used in the 
~^ pK^TT^ reconstruction. 



3 Evaluation of the ratios fs/{fu + fd) and 

fAt/ ifu + fd) 

Perturbative QCD calculations lead us to expect the ratios fs/ (fu+fd) and fA^/{fu+fd) to 
be independent of rj, while a possible dependence upon the b hadron transverse momentum 
Pt is not ruled out, especially for ratios involving baryon species [18]. Thus we determine 
these fractions in different px and r] bins. For simplicity, we use the transverse momentum 
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of the charmed hadron-// pair as the px variable, and do not try to unfold the b hadron 
transverse momentum. 

In order to determine the corrected yields entering the ratio fs/ {fu + fd), we determine 
yields in a matrix of three t] and five bins and divide them by the corresponding 
efficiencies. We then use Eq. |5, with the measured lifetime ratio (r^- + r^o)/2r^o = 
1.07 ± 0.02 yy to derive the ratio fs/{fu + fd) in two rj bins. The measured ratio is 
constant over the whole rj-pT domain. Figure 10 shows the fg/ {fu + fd) fractions in bins 
of pt in two 7] intervals. 
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Figure 9: Efficiencies for D^fi uX, D^fi uX, Dffi uX, A'^fi vX as a function of 7] and 
Pt- 

By fitting a single constant to all the data, we obtain fs/ifu + fd) = 0.134±0.004;^q;q}q 
in the interval 2 < 77 < 5, where the first error is statistical and the second is systematic. 
The latter includes several different sources listed in Table |2} The dominant systematic 
uncertainty is caused by the experimental uncertainty on B{D^ K^K~7t~^) of 4.9%. 
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Figure 10: Ratio between and light B meson production fractions as a function of the 
transverse momentum of the D^fi~ pair in two bins of t]. The errors shown are statistical 
only. 



Table 2: Systematic uncertainties on the relative B^ production fraction. 

Source Error (%) 

Bin-dependent errors 1.0 

B{D^ ^ K-n+) 1.2 

B{D+ 1.5 

bId+ K-K+TT+) 4.9 

semileptonic decay modelling 3.0 

Backgrounds 2.0 

Tracking efficiency 2.0 

Lifetime ratio 1.8 

PID efficiency 1.5 

5° ^ D'>K+Xfi-u +ti 

B{{B-, 5°) ^ DfKXfi-V) 2.0 
^Ibtai 



Adding in the contributions of the and branching fractions we have a systematic 
error of 5.5% due to the charmed hadron branching fractions. The -B^ semileptonic mod- 
elling error is derived by changing the ratio between different hadron species in the final 
state obtained by removing the SU(3) symmetry constrain, and changing the shapes of the 
less well known D** states. The tracking efficiency errors mostly cancel in the ratio since 
we are dealing only with combinations of three or four tracks. The lifetime ratio error 
reflects the present experimental accuracy [1] . We correct both for the bin-dependent PID 
efficiency obtained with the procedure detailed before, accounting for the statistical error 
of the calibration sample, and the overall PID efficiency uncertainty, due to the sensitivity 
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Figure 11: /+//o as a function of for ri={2,3) (a) and ?7=(3,5) (b). The horizontal hne 
shows the average value. The error shown combines statistical and systematic uncertain- 
ties accounting for the detection efficiency and the particle identification efficiency. 



to the event multiplicity. The latter is derived by taking the kaon identification efficiency 
obtained with the method described before, without correcting for the different track mul- 
tiplicities in the calibration and signal samples. This is compared with the results of the 
same procedure performed correcting for the ratio of multiplicities in the two samples. 
The error due to 5° — )■ X n^V is obtained by changing the RS/WS background 

ratio predicted by the simulation within errors, and evaluating the corresponding change 
in fa/ {fu + fd)- Finally, the error due to {B^ , B^) KXn~V reflects the uncertainty 

in the measured branching fraction. 

Isospin symmetry implies the equality of fd and /„, which allows us to compare /+//o = 
^corr(-D"'"//)/ncorr(-D°/i) with its expccted value. It is not possible to decouple the two ratios 
for an independent determination of fu/ fd- Using all the known semileptonic branching 
fractions [Ij, we estimate the expected relative fraction of the and modes from 
5+/° decays to be /+//o = 0.375 ± 0.023, where the error includes a 6% theoretical 
uncertainty associated to the extrapolation of present experimental data needed to account 
for the inclusive h — )■ cji~V semileptonic rate. Our corrected yields correspond to f+/ fo = 
0.373 ± 0.006 (stat) ± 0.007 (eff) ± 0.014, for a total uncertainty of 4.5%. The last 
error accounts for uncertainties in B background modelling, in the D^K^n'V yield, the 
D^pfi'T/ yield, the and branching fractions, and tracking efficiency. The other 
systematic errors mostly cancel in the ratio. Our measurement of /+ / /o is not seen to be 
dependent upon px or rj, as shown in Fig. 11, and is in agreement with expectation. 

We follow the same procedure to derive the fraction fAt/{fu + fd), using Eq. [T] and 
the ratio (r^- + Tjjo)/(2ryio) = 1.14 ± 0.03 [1]. In this case, we observe a px dependence 
in the two rj intervals. Figure 12 shows the data fitted to a straight line 

-Jj^ = a[l + bxp^iGeV)]. (8) 

Ju + Id 

Table [3] summarizes the fit results. A corresponding fit to a constant shows that a 
Pt independent fA^/ifu + fd) is excluded at the level of four standard deviations. The 
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Figure 12: Fragmentation ratio fAt/ifu + fd) dependence upon p'r{A'^fi ). The errors 
shown are statistical only. 



systematic errors reported in Table [3] include only the bin-dependent terms discussed 
above. 

Table |4] summarizes all the sources of absolute scale systematic uncertainties, that in- 
clude several components. Their definitions mirror closely the corresponding uncertainties 
for the fs/ifu + fd) determination, and are assessed with the same procedures. The term 

— )■ D^pXfi~u accounts for the uncertainty in the raw D^pXfi~u yield, and is evalu- 
ated by changing the RS/WS background ratio (1.4±0.2) within the quoted uncertainty. 
In addition, an uncertainty of 2% is associated with the derivation of the semileptonic 
branching fraction ratios from the corresponding lifetimes, labelled Fgi in Table |4} The 
uncertainty is derived assigning conservative errors to the parameters affecting the chro- 
momagnetic operator that influences the B meson total decay widths, but not the A^. By 
far the largest term is the poorly known B{A'^ — )■ pK~7i~^); thus it is quoted separately. 

Table 3: Coefficients of the linear fit describing the pT^A^fi^) dependence of fA^/{fu + fd)- 
The systematic uncertainties included are only those associated with the bin-dependent 
MC and particle identification errors. 



1] range 


a 


b 


2-3 


0.434±0.040±0.025 


-0.036±0.008±0.004 


3-5 


0.397±0.020±0.009 


-0.028±0.006±0.003 


2-5 


0.404±0.017±0.009 


-0.031±0.004±0.003 
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Table 4: Systematic uncertainties on the absolute scale of fh^l{fu + fd)- 



OUUi L-c 




J_>lli U-tJUdiU-CilL CiiUio 


9 9 




9 n 


i\/r x/^i Jii" 

Monte Uarlo modelling 


"1 C\ 

1.0 


Backgrounds 


3.0 


Tracking efficiency 


2.0 


Tsi 


2.0 


Lifetime ratio 


2.6 


PID efficiency 


2.5 


Subtotal 


6.3 


B{A+ pK-7T+) 


26.0 


Total 


26.8 



In view of the observed dependence upon we present our results as 

(Pt) = (0.404±0.017±0.027±0.105) x [1 - (0.031 ±0.004±0.003) x pT(GeV)], 

(9) 

where the scale factor uncertainties are statistical, systematic, and the error on Bi^A^ — ?■ 
pK^n^) respectively. The correlation coefficient between the scale factor and the slope 
parameter in the fit with the full error matrix is —0.63. Previous measurements of this 
fraction have been made at LEP and the Tevatron [3]. LEP obtains 0.110±0.019 [2]. This 
fraction has been calculated by combining direct rate measurements with time-integrated 
mixing probability averaged over an unbiased sample of semi-leptonic b hadron decays. 
CDF measures + fd) = 0.281 ± 0.0121°;°^^^°;^^^, where the last error reflects the 

uncertainty in B{A^ — )• pK^n^). It has been suggested [3] that the difference between the 
Tevatron and LEP results is explained by the different kinematics of the two experiments. 
The average pt of the A'^fi~ system is 10 GeV for CDF, while the 6-jets, at LEP, have 
p 40 GeV. LHCb probes an even lower b pt range, while retaining some sensitivity in 
the CDF kinematic region. These data are consistent with CDF in the kinematic region 
covered by both experiments, and indicate that the baryon fraction is higher in the lower 
Pt region. 
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4 Combined result for the production fraction fs/fd 
from LHCb 



From the study of h hadron semileptonic decays reported above, and assuming isospin 
symmetry, namely fu = fd, we obtain 

^ ) = 0.268 ± 0.008(stat)t[!-[!^^(syst), 

JdJ s\ 

where the first error is statistical and the second is systematic. 

Measurements of this quantity have also been made by LHCb by using hadronic B 
meson decays [1]. The ratio determined using the relative abundances of 5° — )• D^n^ to 
5° ^ D+K- is 



0.250 ± 0.024(stat) ± 0.017(syst) ± 0.017(theor) 



hi 

while that from the relative abundances of -B^ -D^vr" to — t- D^tt~ is 
jA = 0.256 ± 0.014(stat) ± 0.019(syst) ± 0.026(theor). 

JdJ h2 

The first uncertainty is statistical, the second systematic and the third theoretical. The 
theoretical uncertainties in both cases include non-factorizable SU(3)-breaking effects and 
form factor ratio uncertainties. The second ratio is affected by an additional source, 
accounting for the PF-exchange diagram in the B^ — ?■ D^7t~ decay. 

In order to average these results, we consider the correlations between different sources 
of systematic uncertainties, as shown in Table [5j We then utilise a generator of pseudo- 
experiments, where each independent source of uncertainty is generated as a random 
variable with Gaussian distribution, except for the component 5° — )■ D^K^ ^^u^X , which 
is modeled with a bifurcated Gaussian with standard deviations equal to the positive and 
negative errors shown in Table [5j This approach to the averaging procedure is motivated 
by the goal of proper treatment of asymmetric errors [21] . We assume that the theoretical 
errors have a Gaussian distribution. 

We define the average fraction as 

fs/fd = ai(/.//d)si + a2(/.//d)hi + a3(/s//d)h2, (10) 

where 

«! + 02 + «3 = 1- (11) 

The RMS value of fs/fd is then evaluated as a function of ai and 02- 

We derive the most probable value fs/fd by determining the coefficients ai at which 
the RMS is minimum, and the total errors by computing the boundaries defining the 
68% CL, scanning from top to bottom along the axes ai and a2 in the range comprised 
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Table 5: Summary of the systematic and theoretical uncertainties in the three LHCb 
measurements of fsl fd- 



Source 




Error (%) 








(/s//d)sl 


ifs/ /d)hl 


{fs/ fd)h2 




Bin dependent error 


1.0 






Uncorrelated 


Semileptonic decay modelling 


3.0 


_ 


_ 


Uncorrelated 


Backgrounds 


2.0 


_ 


_ 


Uncorrelated 


Fit model 




2.8 


2.8 


Uncorrelated 


Trigger simulation 




2.0 


2.0 


Uncorrelated 


Tracking efficiency 


2.0 






Uncorrelated 




+4.1 
-1.1 






Uncorrelated 




2.0 






Uncorrelated 


Particle identification calibration 


1.5 


1.0 


2.5 


Correlated 


B lifetimes 


1.5 


1.5 


1.5 


Correlated 


B{D+ K+K-TT+) 


4.9 


4.9 


4.9 


Correlated 


B{D+ K-TT+n-) 


1.5 


1.5 


1.5 


Correlated 


SU(3) and form factors 




6.1 


6.1 


Correlated 


l^-exchange 






7.8 


Uncorrelated 



between and 1. The optimal weights determined with this procedure are ai = 0.73, and 
a2 = 0.14, corresponding to the most probable value 

= 0.267in^S. 

The most probable value differs slightly from a simple weighted average of the three 
measurements because of the asymmetry of the error distribution in the semileptonic 
determination. By switching off different components we can assess the contribution of 
each source of uncertainty. Table [6] summarizes the results. 



5 Conclusions 

We measure the ratio of the B^ production fraction to the sum of those for B~ and B^ 
mesons fs/ (fu+fd) = 0.134±0.0041o;qJq, and find it consistent with being independent of r] 
and pt- Our results are more precise than, and in agreement with, previous measurements 
in different kinematic regions. We combine the LHCb measurements of the ratio of B^ to 
production fractions obtained using b hadron semileptonic decays, and two different 
ratios of branching fraction of exclusive hadronic decays to derive fs/fd = 0.2671q;q2o- 
The ratio of the A° baryon production fraction to the sum of those for B~ and B^ mesons 
varies with the px of the charmed hadron muon pair. Assuming a linear dependence up 
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Table 6: Uncertainties in the combined value of fs/fd- 



Source 


Error (%) 


Statistical 


2.8 


Experimental systematic (symmetric) 


3.3 


^3(5° ^ D^K+Xfi-u) 


+3.0 
-0.8 


B{D+ K~7T+7r~) 


2.2 


BiDf -> K+K-7r+) 


4.9 


B lifetimes 


1.5 


B{ByB- DfKXfi-u) 


1.5 


Theory 


1.9 



to pt = 14 GeV, we obtain 

= (0.404 d= 0.017 d= 0.027 ±0.105) x [1 - (0.031 ± 0.004 ±0.003) x pT(GeV)], (12) 

Ju + Jd 

where the errors on the absolute scale arc statistical, systematic and error on B{A'^ — >■ 
pK~Tr~^) respectively. No r] dependence is found. 
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